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ABSTRACT: 



The numerical frontogenesis model of Williams (1972) is modified 
to include horizontal and vertical turbulent diffusions of heat and 
momentum* The turbulent diffusions are represented with constant 
coefficients, and an Ekman layer is added to the basic deformation 
field. The numerical solutions show realistic quasi-steady fronts 
forming within 1-2 days. These solutions are examined and 
compared over a wide range of the various coefficients. 
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1. Introduction 



Hoskins and Bretherton (1972) have shown analytically, and Williams 
(1967, 1972) has shown numerically that discontinuous fronts can form within 
a finite period of time if no turbulent diffusion is present. These 
studies suggest that a discontinuity will form within 24 to 36 hours when 
reasonable initial conditions are used. If turbulent diffusion is present, 
it can be expected that a balance will be achieved between the fronto- 
genetic advections and the turbulent diffusions of heat and momentum. 

The front should remain in this state of quasi-balance as long as the 
large scale deformation field causes frontogenetic advections around the 
front. The purpose of this study is to obtain and examine steady-state 
frontal solutions. Welander (1963) obtained steady-state frontal solutions 
with a boundary layer technique, but his formulation did not involve 
a frontogenetic large-scale deformation field. 

The turbulent diffusion coefficients for heat and momentum are taken 
to be constant in this study, although these coefficients are actually 
functions of the wind and temperature fields. Smagorinsky (1963) and 
Leith (1968) have proposed forms for the horizontal diffusion coefficients 
which are dependent on the local wind field. Many studies have treated the 
form of the vertical diffusion coefficients near the surface (Estoque 1960, 
Deardorff 1970, 1972). In this study we will carry out a large number of 
integrations with different values of the various constant coefficients. 

The range of values will include the mean values which are given by the 
various nonlinear formulations. 

This investigation will employ the numerical model of Williams (1972); 
this reference will hereafter be referred to as W72. The physical model is 
essentially the same as one of the models treated by Hoskins (1971) and 
Hoskins and Bretherton (1972). The frontogenesis is forced by a non- 
divergent horizontal wind field which contains stretching deformation. 

The hydrostatic primitive equations are used with the Boussinesq approximation. 
In the model the time dependent quantities are functions of y and z only. 
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The basic deformation wind field is constant in time and is independent 
of height, except in the surface boundary layer. 

In Section 2, the basic forecast equations are developed, and a 
simplification which keeps the problem two-dimensional is discussed. 

The initial conditions are presented in Section 3. The solutions which 
contain only horizontal diffusion are discussed in Section 4, and those 
which contain only vertical diffusion are discussed in Section 5. The com- 
bined experiments which include both types of diffusion are examined in 
Section 6 and the conclusions are given in Section 7. 



- 2 - 



2 . Basic equations 



The hydrostatic Boussinesq eqtiation with diffusion niay be written: 



at 

at 



+ V 



+ V . 



OK) + (wV) + V(|) + Qc X V = Av2 V + 

(0V) + ^ (w0) = AV^0 + Cg a^0 + Qa 
^2 ^ 



c a^v 

■m 




( 2 . 1 ) 

( 2 . 2 ) 



V • K + = 0, 

^ = ge ^ 
az 0Q 



(2.3) 

(2.4) 



where the notation is the same as used in W72. The quantity A is the 

horizontal turbulent diffusion coefficient; note that the same value is 

used for both heat and momentum. The quantities C and are the vertical 

me 

turbulent diffusion coefficients for momentum and heat respectively. 

A convective adjustment process is represented by the function Q • The 

a 

boundary conditions at the top and bottom of the domain are: 
w(x,y,o,t) = w(x,y,H,t) = 0 
V(x,y,o,t) = 0 

Sn ^ (x,y,H,t) = 0 

bz 

(x,y,o,t) = ^0 ^ (x,y,H,t) = 0 

az Sz 

The lower boundary conditions are realistic for fronts over land, but 
the upper lid boundary conditions give only a rough approximation to the 
tropopause. Upper frontogenesis has been treated by Hoskins and Bretherton 
(1972) and Hoskins (1972). In this paper our main concern is with fronts 
in the lower troposphere. 
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The following fields which satisfy the boundary conditions (2.5) 
are approximate solutions to the steady~state equations [ (2.1) ~ (2.4)] . 



( 2 . 6 ) 



( -Bz “bz 

V = U = D I [x(l -e cos Bz) + ye sin BzJ4> 

■[y(l - e cos Bz)-xe sin Bz] J | 
w = 0 

^ = -D^(x^ + y^)/2 - fDxy 

9=0 

where B = (f/2Cni)^ where D and f are constant. These relations 
satisfy the steady-state equations if some small advection terms in the 
boundary layer are neglected. These small terms are usually dropped in the 
Ekman theory and they are of the order of the Rossby Number which for 
these fields is D/f. In this paper o/f = 0.1 for all cases. 

We now subdivide our dependent variables as follows: 

V = U (x,y,z) + u(y, 2 , t)^ + v (y,z,t) ^ 
w = w(y, z, t) 

0 = 0(y, z, t) 
f = $(x,y) + rr(y, z, t) 

It is noted that all departures from the steady deformation solutions (2.6) 
are assumed to be independent of x. If we substitute the relations (2.7) 
in the ;L component of (2.1) we obtain: 



(2.7) 



^ ^ (uv) + ^ (wu) H- uD(l - 

St Sy Sz 

-D[y(l - e cos Bz ) - xe sin Bz] ^ + vD e sin Bz 

-Bz 



4wD^ r xe"'"'' (1- cos Bz) + ye'“^ sin Bz 

5z L 



-Bz 



cos Bz) 

by 



( 2 . 8 ) 






2 2 
fv + kh u + C mo u 



Please note that the assumption that u is independent of x for all time 
is violated in (2.8) in the two terms where x appears. This indicates 
that the other quantities which were assumed independent of x are, in fact. 
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x-dependent and this would bring other terms into the equations. In this 
study we will neglect these added terms and apply (2.8) and the other 
equations at x = 0. If the initial u, v, w, 0, and rr fields are independent 
of X the error from this approximation will grow slowly and it will be 
confined to the boundary layer since the x-dependent terms in (2.8) are 
zero outside the boundary layer. Since frontogenesis occurs very rapidly, 
it is expected that the development of x-variations in the dependent 
variables would have only a small effect on the resulting quasi-steady 
front. In any case these effects could not be observed in the atmosphere, 
because atmospheric fronts always have some variation in the basic fields 
along the front which would be much more important. 

When Eq. (2.8) is evaluated at x = 0 it becomes: 



5u + S(uv) + 5(wu) - Fu ^ + rv5u - vSv e sin Bz 
St Sy Sz Sy Sy Sy 

-VwS ( e sin B^=fv + AS^u + C mS^u , 

Sz ^^2 ( 

where; V = -Dy and r = 1 - e'^^cos Bz. The j component of (2.1) applied 
at X = 0 takes the form: 



— B z • 

Sv + ^(vv) + ^ (wv) + (vV) - uSv e sin Bz 
St Sy Sz Sy Sy 

-wV5 (e ^^c©s Bz )= - fu + . (2.10) 

^ ^ Sy2 ^ 



In a similar manner Eq. (2.2) may be written: 

^ + ^(v9) + ^(W0) + ryM = + 00^20 + , (2.11) 

St Sy Sz Sy Sy2 Sz2 

and Eqs. (2.3) and (2.4) become: 

= 0 , ( 2 . 12 ) 

(2.13) 



Sv + 



Sz 



= 

06 ^ 
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The boundary conditions at the top and the bottom of the domain are 
given by: 

u = V = w = 00^9 = 0 , z = 0 

5z , (2.14) 

Cm^ = Cm^ = w = CqS 0 = 0 , z = H 
3z 5z 5z 

where H is the distance between the rigid plates. The no slip condition at 
the lower boundary and the no strex condition at the top roughly represent 
the troposphere. 

If we define the vertical average of a quantity as 
H 

<( )> = 1 \ ( )dz, 

H o 

and integrate the hydrostatic equation (2.13) with respect to z and remove 
the vertical mean, we obtain 

z z 

0 dz - < ^ 0 dz >"| . (2.15) 

o 

Take the vertical average of (2.12) and use the boundary conditions (2.14), 
which gives 

5 <v> = 0 . 

Sy 

This equation states that the total mass flux in the j direction is inde- 
pendent of y. In W72 a symmetry argument was used to show that this flux 
must be zero. This argument does not hold strictly in our case, but it 
can be expected that at large distances from the frontal zone the distur- 
bance mass flux will vanish. 

Thus we set 

<v> = 0 . (2.16) 
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If we utilize the development in W?2 we can rewrite (2.10) in the following 
form: 

5v + ^ [vv “ <w>] + 5(wv) + ^ (pW - <pv>V) 
bt ^ bz by 

- ^ ^^sin Bz -<ue sin Bz>] - v|w5 (e cos Bz) 

by L Sz 

~<wb Bz)>1 = - 5 (n - <n>) - f(u - <u>) 

Bz J ^ 

+ AB^v + Cm r B^v + 1^ \ "I 

^ hVSz^^oJ' ^2.17) 

The Eqs. (2.9), (2.11), (2.12), (2.15), and (2.17) forma complete 
set which can be solved by a pure marching process. The finite difference 
equations conserve mean squares in the advection terms and are described 
by Williams (1967). The frictional terms are evaluated at the previous 
time step to avoid computational instability. Computational boundaries 
are introduced at y = ± Y and these are the same as used in W72. On these 
boundaries u and 0 are held fixed and v = 0. 

3. Initial conditions 

The initial conditions are the same as those employed in W72. The 
initial temperature field is given by 

0(y, z, o) = (z - h/2) - a(2/iT) arctan (sinh ay), (3.1) 

Bz 

where a = f tt H (g0o^ ^^/Bz) . The quantity B0j/Bz, which is constant 

is the initial static stability and a is one half the total horizontal 
temperature variation 

The initial x - component of the velocity is given by 

u(y, z, o) = 2 gag (z - H/2) sech ay. (3.2) 

n f0Q 

This equation is obtained by substituting (3.1) into the thermal wind 
equation and integrating. In this study the v field is the divergent part 
part of the wind. We initialize the v field by equating it to the diver- 
gent part of the wind which is obtained from the quasi-geostrophic equations 
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The details are given in W72* 

It is noted that these initial conditions do not satisfy the boundary 
conditions on u, v, and 0 at the upper and lower boundaries. A period of 
adjustment will be required to form the surface friction layer and some 
oscillations may be observed later. 

4. Numerical experiments with horizontal diffusion 

All numerical results to be shown use the following values for the 
constants : 

f = 10"^ sec"^ , £ = ,0327 m sec"2(OK)‘l 

6o 

H = 9 km, Y = 1800 km 
D = 10 ^sec ^ km“^ 

a = 12.56K 

These values were also used in W72. 

In this section we examine the numerical solutions which include only 
horizontal diffusion, so that Cm = Cq = 0. Many mamerical experiments 
are examined in this section which have values of the horizontal diffusion 
coefficient ranging from A = 10“^^ m^sec*^ to A = 10^m^sec“^. All experiments 
in this section have a vertical increment of Az = 333 m. The horizontal 
mesh length is either 20 km or 60 km depending on the expected frontal 
scale. The convective adjustment process Qa is not required in these 
experiments. 

A reasonable measure of the width of the frontal zone is given by 
d = e(-Y. Z. o) - e(Y. z, o) , 

(4.2) 

max 

where 50 /5y is approximated by a one-sided difference. Figure 1 shows the 
variation of d for 3 different experiments at the lowest level (z = 167 m) . 
The lowest curve which is taken from W72 is for the frictionless case 
(A = 0) with Zy = 20 km. In this experiment no steady-state is reached 
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although d is bounded by [sy because of the truncation error. In W72 it 
was estimated that a discontinuity would occur at t = 1.75 days if there 
were no numerical truncation error. The middle curve corresponds to the 
experiment with A = 3 x 10^m^sec“^ and A7 = 20 km, and shows that the 
solution reaches a steady-state at 2 days and is close to the limiting 
scale within 1.5 days. The upper curve is for A = 5 x 10^ m^sec*"^ and 
Ay = 60 km. This experiment shows the same behavior as with the smaller 
diffusion coefficient except that the limiting scale is much larger. 

Figure 2 contains cross sections of 0^ = 0 - , where 0 is the 

horizontally averaged initial potential temperature. Only the lower 
half plane is shown because the solutions of this section satisfy the sym- 
metry condition which is given by Eq. (5.4) in W72. Also, the outer 
portion of the domain between lyj.= 1200 km and jy] = 1800 km is not 
shown. The solutions shown in the remainder of this section are for 
t = 4 days. The solutions in Fig. 2 may be compared with the initial 
conditions and the frictional solutions which are shown in Figure 4 
of W72. The solution for A = 3 x lO^m^sec"^ is very similar to the 
frictionless solution at t = 1.25 days except that the frictional case 
has a larger tilt. The solution for A = 5 x lO^m^sec"^ is similar to the 
solution for A = 3 x 10^m2 sec^ except that the gradients are much 
weaker and the frontal zone has slightly less tilt. 

Figure 3 contains cross sections of u for the same cases and these 
can be compared with Figure 5 in W72. As noted above, the solution for 
A = 3 X 10^m2 sec”^ is very similar to the frictionless solution at t = 1.25 
days, except that the tilt is greater for the frictional case. The solution 
for A = 5 X lO^m^sec"^ is similar to the solution A = 3 x lO^m^sec"^ 
except that the velocities are much smaller and the gradients are much 
weaker. 

Figure 4 summarizes all the numerical experiments which contain 
only horizontal diffusion. The figure contains the frontal scale at 
z = 167 m as a function of the diffusion coefficient on a log-log plot. 

The dots represent experiments with Ay = 60 km and the circles represent 
experiments with Ay = 20km. For comparison, let us derive the quasi- 
geos trophic solution at z = 0. The steady-state first law of thermodynamics 
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at the surface is 



-Dy^'= 

Sy ^ * 



The solution to this equation is 
6'(y, o) = |e'(“,o) - 0'(-oo,o)| (2 Att/d) ^ 

O ^ ' 

When this form is substituted into the definition we obtain 



^ exp/-D^j 



dy 



d = (2 ATt/d)^® , 



(4.3) 



(4.4) 



(4.5) 



which is shown in Figure 4. The plotted data approach the line (4.5) 
as A becomes large, which is expected because the larger scales involved 
should lead to more nearly quasi-geostrophic solutions (Williams, 1967). 

As A is decreased the points fall ferliier below the curve although the 
slope is not much greater than )/ . This is the effect of the advections 
of temperature and momentum by the divergent part of the wind. These 
effects are neglected in the quasi-geostrophic equations. For the smallest 
values of A the plotted data again approach and cross the quasi-geostrophic 
curve. This is a consequence of the truncation error which arises when 
the scale becomes very small. 

5 . Numerical experiments with vertical diffusion 

In this section we examine the numerical solutions which include only 

the vertical diffusion so that A = 0. Both vertical diffusion coefficients 

2-1 2-1 

are varied over a range from 1.25 m sec to 15 m sec . The grid increments 
for all of these experiments are Ay = 60km and Az = 200m. This value of 
the vertical increment was chosen to resolve better the vertical structure 
of the frontal zone and the surface friction layer. In these experiments 
the convective adjustment process is required at times near the surface 
on the cold air side of the front. 

In most of the experiments and Cq are not equal because they are 
normally different in the lower troposphere due to the influence of thermal 
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stratification. Figure 5 shows the variation of the frontal scale for 

3 combinations of and Cq. These curves may be compared with Figure 1 

which gives solutions with horizontal diffusion only. The solutions in 

Figure 5 all show oscillations which were not present in the horizontal 

diffusion cases. These oscillations arise from the development of the 

surface boundary layer which is not present in the initial conditions. 

The interior fields must adjust to the no-slip velocity condition and to the 

insulated potential temperature condition at the lower boundary. The 

delay in frontogenesis in the first few hours between Figures 5 and 1 

is also due to the adjustment process. At t = 4 days the upper 2 curves 

in Figure 5 are nearly steady while the lowest curve (Cm = 10m2sec”^, 

2 -1 

Cg = 1.25m sec ) still contains oscillations. We will compare the 
various solutions in this section at t = 4 days although those experiments 
which have a small Cg may not be quite steady. 

Figure 6 contains cross sections of 0^ for the 3 experiments at 
t = 4 days. The solutions in this section do not have the vertical symmetry 
properties which were noted in the previous section because the upper and lower 
boundary conditions are not equal. Nevertheless we display only the 
solutions for the lower half of the domain because we are concerned mainly 
with fronts near the surface and also because the solutions near the upper 
boundary are not very realistic. The solutions in Figure 6 may be compared 
with the horizontal diffusion solutions which are given in Figure 2. 

Comparison of these figures shows some important differences between 
steady-state fronts in which only vertical diffusion is present. The 
fronts with vertical diffusion all show a greater tilt than those with only 
horizontal diffusion. The solution with the smallest vertical heat dif- 
fusion coefficient (Cm = 10 m^sec"^, Cg = 1.25 m^sec"^) has a large tilt 
in the lower 2 kilometers of the atmosphere. The vertical diffusion in 
these experiments causes the large gradient region to extend higher into 
the atmosphere than for those experiments with only horizontal diffusion. 

It is noted that Sfl '/Sz is negative near the lower boundary, and this comes 
from the insulated boundary condition at z = 0: 

^ = 0 = ^1 + ^' 

Sz 5z Sz ’ 

since S^j/Sz is positive. 
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Figure 7 contains cross sections of the u field for the 3 experiments 
at t = 4 days. These solutions can be compared with the horizontal dif- 
fusion solutions which are shown in Figure 3. All of the solutions show 
the presence of the surface friction layer with the no-slip boundary 
condition. The thickness of the surface boundary layer appears to vary 
in proportion to CniVs as would be predicted by Ekman theory. When these 
fields are compared with the temperature fields in Figure 6, it can be 
seen that the regions of large cyclonic shear and large vertical shear 
occur in the zones of large temperature gradient. As would be expected, 

the maximum speed occurs in the experiment with the minimum value of the 

2-1 2-1 

momentum diffusion coefficient (Cm = 1.25 m sec , Cg = 5 m sec ). 

Figure 8 summarizes all the numerical experiments which contain only 

vertical diffusion. The data for each experiment at t = 4 days and z = 100 m 

are plotted as a function of Cm and Cg. The number above each point is 

the frontal scaled and the number below each point is the frontal location 

y^. A striking feature of the diagram is the tendency for d to remain 

relatively constant as is increased when Cg is held fixed. Conversely 

d increases rapidly with Cg when Cm is held fixed. The quantity y^ is 

a measure of the mean frontal tilt because in the center of the domain 

(z = 4.5 km) the frontal zone is located near y = 0 (see Figures 2 and 6). 

As noted above some experiments have a much larger tilt in the lower 2 

kilometers. Figure 8 shows a slow increase in the tilt as Cm is increased 

for a fixed Cg. When C^ is held fixed the tilt decreases more rapidly 

as Cg is increased. The experiments for the special cases C = 0 and 

2 -1 

Cg = 0 are not included in the figure. The experiment for = 10 m sec 
and Cp = 0 contains large special oscillations while the solution for Cm = 0 
and Cg = 5 m sec contains large time fluctuations. 

Some of the gross features of Figure 8 may be interpreted as follows. 
Consider the behavior for fixed Cg. As Cj^^ is increased we expect the 
vertical wind to decrease and through the thermal wind equation this should 
decrease the horizontal temperature gradient. There is, however, an 
opposing effect in the surface boundary layer. In the region of cyclonic 
vorticity we expect horizontal mass convergence in the boundary layer; 
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this should occur even if the boundary layer is not of the pure Ekman type 
(Mak, 1972). This convergence will tend to build up the temperature 
gradients within the friction layer. This effect is opposed in part by 
the vertical diffusion of heat which is proportioned to Cq. As is 
increased the boundary layer will become thicker, as was pointed out earlier 
in this section. This increase in the thickness of the friction layer 
will tend to increase the temperature gradients as long as is held 
fixed. The numerical solutions show that this frictional convergence 
effect on the temperature gradient is approximately balanced by the momen- 
tum diffusion effect through the thermal wind. Therefore d remains 
nearly constant as is increased with Cg held fixed. 

When C is held fixed and Cn is decreased the vertical scale of the 
m o 

frontal zone becomes very small. This can be seen in Figure 6 if the lower 
2 cross sections are compared. The horizontal temperature gradient is 
somewhat constrained through the thermal wind relation by vertical momentum 
diffusions. The decrease in the vertical scale, which comes from a de- 
crease in Cg, must be associated with an increase in the tilt of the 
frontal zone if the change in the horizontal gradient is to be minimized. 

This behavior is clearly present in Figure 8. 

6. Numerical experiments with both horizontal and vertical diffusion 

In this section we will examine the solutions which contain both 
horizontal and vertical diffusion. Individual solutions will not be 
shown because they are essentially combinations of the solutions shown in 
earlier sections. The principal features of the solutions are summarized 
in Tables 1 and 2. All solutions have space increments of Zy = 60 km 
and = 200 m. The data are extracted at t = 4 days, and all of the 
data in the tables are for the lowest level (z = 100 m) . 

Table 1 contains the frontal scale d for various combinations of the dif- 
fusion coefficients. The first column gives the values for sero vertical 
diffusion and the first row gives the values for zero horizontal diffusion. 

The smallest value of d occurs at the lowest level for each experiment 
except for the two with footnotes. In the two exceptional cases, the 
minimum occurs at the second point and this is caused by the vertical 
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diffusion of momentum. It may be seen that the value of d for each 
experiment is a little less than the sum of the values obtained from the 
experiments where one of the coefficients is zero. 

Table 2 contains the frontal locations for various combinations of 
the diffusion coefficients. In the first column, which is for zero ver- 
tical diffusion, the tilt decreases slightly with increasing horizontal 
diffusion coefficient. The combined experiments all show values which 
are closer to the A = 0 than to the Cjjj = Cq = 0 solutions. However, 
for each combination of and Cg, the tilt does decrease as A is in- 
creased. We conclude that the frontal scale is determined by both dif- 
fusion processes, but that the frontal tilt is determined principally 
by the vertical diffusion. 

7. Conclusions 

The numerical model of Williams (1972) is modified in this study to 
include horizontal and vertical diffusions of momentum and heat. The 
turbulent diffusions are represented with constant coefficients. An Ekman 
layer is added to the basic deformation field and some other terms which 
are related to the modified deformation field. The numerical solutions 
become quasi-steady within 1-2 days. These results suggest that 
atmospheric fronts form in a fairly short period of time (1-2 days) 
and then move in a state of quasi-balance where the frontogenetic processes 
are balanced by diffusion processes. The time scale of this second phase 
would be related to the time scale of the deformation field which caused 
the frontogenesis. 

The solutions with only horizontal diffusion are similar to the transient 
solutions obtained in W72 at particular times. The horizontal scale at 
the lowest level follows approximately an A^ dependence. The values obtained 
are less than the exact quasi-geostrophic solutions at z = 0 and they approach 
the latter as A increases. For reasonable values of A the scale is close 
to the quasi-geostrophic value, but the structure is very different since 
the quasi-geostrophic solutions contain no vertical tilt. 

The solutions with only vertical diffusion are strongly dependent 
on the ratio of to Cg. All of these solutions show a deeper region 
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of large gradient than for the solutions with only horizontal diffusion. 

When Cg is held fixed, the scale d shows little change as is increased. 

On the other hand, when Cm is held fixed, d increases rapidly with increasing 
Cg. The frontal tilt increases when Cm is held fixed and Cg is decreased. 

The solutions for C^ = 0 or Cg = 0 are not acceptable steady solutions. 

The solutions which include both horizontal and vertical diffusion 
are combinations of the separate solutions. The values for d are a little 
less than the sum of the values from each of the diffusion processes. 

The tilt of the fronts is principally determined by the vertical diffusion 
effects. 

The steady-state frontal solutions obtained in this paper are reasonable 
with respect to observed fronts. The vertical diffusion terms are clearely 
required to give the proper structure near the surface. Both Cm and Cq 

must be included although a value C that is somewhat less than Cm may be 

0 

appropriate. The importance of the horizontal diffusion processes can 
probably be determined only from turbulence studies in frontal zones. 

The solutions found in this paper would be affected if nonlinear 
diffusion coefficients were used, but we expect the same general type of 
solution. However, if C^ were allowed to decrease above the friction layer, 
the maximum gradient might occur higher in the atmosphere. 

The tilt of the solutions obtained would be modified if the whole 
system were allowed to have a net translation. Then it would be possible 
to distinguish between cold and warm fronts. The addition of latent 
heat should be important in determining the structure of many fronts. 
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Table 1. Values of d in kilometers 



Cg 

A 


0,0 


10, 1.25 


10, 10 


1.25, 5 


0 




172 


399 


241 


3 X 10^ 


94 


239 


450 


291 


10^ 


181 


321 


523 


370 


5 X 10^ 


492 


666^ 


2 

828 


644 



1 Minimum at z = 300 m of 631. 

2 Minimum at z = 300 m of 827. 



Table 2. Values of in kilometers 



A 


0,0 


10, 1.25 


10, 10 


1.25, 5 


0 


_ _ _ _ 


-1010 


-820 


-720 


4 

3 X 10 


-545 


-995 


-805 


-715 


10^ 


-540 


-985 


-755 


-700 


5 X 10^ 


-510 


-985 


-715 


-640 
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FIGURE CAPTIONS 



Figure 1. 

Figure 2. 

Figure 3. 

Figure 4, 

Figure 5. 

Figure 6. 
Figure 7. 
Figure 8. 



Time variation of the frontal scale at z = 167 m for experiments with 
only horizontal diffusion. 

5 2 "1 

Cross sections of 0^ = 0 - 0^ for A = 5 x 10 m sec and 
A = 3 X 10^ m^sec*"^ at t = 4 days. 

5 2 -1 4 2 -1 

Cross section of u for A = 5 x 10 m sec and A = 3 x 10 m sec 

at t = 3 days. 

The variation of the frontal scale at z = 167 m as a function of the 
horizontal diffusion coefficient. 

Time variation of the frontal scale at z = 100 m for experiments 
with only vertical diffusion. 

Cross section of 0^ for 3 combinations of C and C at t = 4 days. 

m 0 ^ 

Cross sections of u for 3 combinations of C and C^ at t = 4 days. 

m 0 ^ 

The frontal scale and frontal location as functions of C and 
Cg. These data are for z = 100 m and t = 4 days. ^ 
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Fig. 1 
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Fig. 3 
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Fig. 7 
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Fig. 8 
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